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Abstract 



'^ \ A model is constructed in which pair potentials are combined with the cluster expansion method 

in order to better describe the energetics of structurally relaxed substitutional alloys. The effect of 
structural relaxations away from the ideal crystal positions, and the effect of ordering is described 
by interatomic-distance dependent pair potentials, while more subtle configurational aspects asso- 



i-^ ' dated with correlations of three- and more sites are described purely within the cluster expansion 

O ' formalism. Implementation of such a hybrid expansion in the context of the cluster variation 



method or Monte Carlo method gives improved ability to model phase stability in alloys from 



^ \ first-principles. 
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The lattice gas model has been very effective for modeling substitutional and interstitial 
alloys and compounds.- Although it is a generalization of the classical Ising model, when 
combined with effective interactions extracted from ab initio total energies through the 
cluster expansion method (CEM),^ it provides the basic framework for the modern theory 
of alloys .-'^'^ 

Experience has learned that on perfect lattices the CEM converges rapidly, requiring only 
clusters with a few sites for the thermodynamic modeling of alloy phase stability, see e.g. 
Ref. [Sj. However, when structural relaxations play an important role, as in alloys involving 
constituents with large size differences, convergence of the CEM becomes poor and typically 
long-ranged effective pair and many-body interactions are necessary.- In some alloy systems 
relaxations effects are dominant J In previous calculations this problem was treated by fully 
relaxing structures when performing first-principles calculations and performing the cluster 
expansion not over the internal energy but over other expedient thermodynamic potentials 
such as the enthalpy.- However, this method fails in a number of instances: (i) When the 
terminal phases have different crystal structures, it frequently happens that relaxation of 
the unstable structures leads to the stable structure. For example in the case of bcc and fee 
structures, while intermediate relaxed structures may exist in clearly fee or bcc derived form, 
for the pure endpoints it has been found that when either fee or bcc is stable, unrestricted 
relaxation of the unstable structure is not possible without arriving at the stable structure, 
(ii) Another difficulty is that as a function of temperature, quite apart from configurational 
changes, there are also changes in the relaxed structure as a result of the lattice vibrations.- 
This aspect of the temperature dependence of the relaxation energy is not presented in the 
current implementation of the CEM. (iii) When structures relax it can become impossible 
to uniquely associate a relaxed cluster with a cluster in the unrelaxed structure. For exam- 
ple, an fcc-based ordered structure such as LIq that relaxes to the bcc-based B2 structure 
presents topological difficulties in that it becomes very difficult to properly categorize and 
count the nearest- and second- nearest neighbor pairs. Although for this specific problem an 
intermediate body centered tetragonal structure can be devised that describes both fee and 
bcc crystal structures,i2 generally, when several modes of distortion exists the definition of 
intermediate or generalized underlying crystal structures becomes impractical. These limi- 
tations of the CEM come about because it was originally conceived for fixed perfect lattices 
and not for relaxed structures.^ Current implementations of the CEM give errors that typi- 



cally are of the order of 300 K for order-disorder transition temperatures Tc. Kikuchi claimed 
that lattice distortion is the main reason for the discrepancy between ab initio-hased and 
experimental Tc.-^ As an illustration one can mention the Tc of LI2 NisAl: when unrelaxed 
structures are used in the CEM, a fortuitous agreement with experiment is found.— »i^ When 
relaxed structures are used in CEM the agreement with experiment worsens^^ and only a 
setter description of structural relaxations can improve agreement, see e.g. Refs. |9| and 
14l |. For simulations of kinetics also, relaxation from the ideal lattice sites may greatly 
influence the energy barriers for atom-exchanges. Therefore, current kinetic Monte Carlo 
diffusion modeling with local-environment independent parameters might be made more re- 
alistic by explicitly treating relaxations. A proper accounting for local distortions might 
much improve the change of the melting temperature under high pressure also.— 

Here, it is our aim to treat the effects of relaxations by proposing an efficient energy func- 
tional. Firstly, an energy function that depends on atomic position must be constructed, 
which is completely different from the usual CEM scheme. Obviously, a pair potential model 
provides the simplest and most widely used function with this characteristic. Although pair 
contributions are the largest energy contribution in alloys and compounds, effective multi- 
body interactions are necessary for accurately describing cluster occupation competitions.-*^ 
However, atomic position dependent multi-body potentials are difficult to obtain and imple- 
ment, so that we opt to include these effects through the effective cluster interactions (ECI) 
as efficiently generated by the CEM. Clearly, it is attractive to combine the CEM and the 
pair potential approach to model the lattice distortion energy. 

Assuming that the magnitude of the distortions is moderate, we may Taylor expand the 
energy as a function of the atomic positions in the vicinity of the relaxed atomic coordinates 
up to second order. This means that we assume that the relaxation energy is mainly due to 
the pair contributions.->i^ This is visualized in Figure [1] where each atom is allowed to move 
(continuously or discretely) within a volume Q around its ideal position.— Thus, in a lattice 
gas treatment each site is characterized not only by its occupation variable but also by its 
displacement vector. Then the energy of the system per atom can be written as 

E =^ Va^a + ^ XI // ^"^"(^«i +^i- Uj)pmn{Ui, Uj)dUidUj, (l) 

a^ij rnn,ij 

where a represents a non-pair cluster, Uj is the displacement vector of atom i from its ideal 
site, Rjj is the vector between unrelaxed sites i and j, and A^ is the number of atoms. The 



FIG. 1: Distortion model for the 2-dimensional square lattice which permits atoms to move around 
their ideal sites within a volume (area) fi centered around the lattice sites (bold dots), as shown 
by the area sprinkled with thin discrete points. Each atom can move within this area without 
constraints derived from symmetry, as exemplified by the arrow. 



summations of m, n in the second term run over atomic species A and B, while ij runs 
over all pairs. The density matrix pm„(uj, Uj) indicates the probability for an "mn" type ij 
atomic pair with site i (j) having a displacement Ui(uj). As for the determination of the 
pair potentials Vmn and the non-pair ECI Va-, we propose two practical schemes: the first is 
based on the lattice inverse method due to Chert^"— and the second employs empirical pair 
potential models. 

Scheme I: A completely ah initio method for finding the pair potential is the lattice 
inverse method based on the modified Mobius inverse^*^ from number theory. It has been 
shown to give reasonable energy functions for a wide range of materials.— First, the cohesive 
energies E{a) for the pure elements A and B with a simple crystal structure /3, such as fee 
or bcc, are computed as a function of the lattice parameter a over a wide range of a values 
with first-principles methods. The pair potential between like atoms, say the AA pair, is 
then extracted from E[a) through an exact transformation^"^ 



vU^) = 2Y,l'{n)Ei[h^{n)x] 



(2) 



71=1 



where the inverse coefficient I^ is derived elsewhere^^ii^s^ and its values for fee and bcc have 
been conveniently tabulated.^® The symbol x in Eq.(l2]) is the nearest neighbor distance, 
while h^{n) are coefficients related to the n^^ coordination shell in the /3 structure, also 
tabulated elsewhere.—"^ Once the pair potentials Vaa and Vbb have been obtained, the pair 
potential related to unlike pairs is determined from the cohesive energies of some unrelaxed 



ordered structures with intermediate composition.— The cohesive energy due to AA and BB 
pairs in the ordered structures as computed with the pair potentials derived for the pure 
elements, is subtracted from the cohesive energy of the ordered structure. This remainder 
of the cohesive energy of the ordered structure is due to the unlike bonds and it is inverted 
just as was done for the pure elements, but of course, now only AB pairs are considered in 
the structure.— 

In practice, there may be more than one elemental structure p. Then, slightly different 
pair potentials may be generated for the different structures. Experience shows the potentials 
derived from different structures to be rather similar,—!^ so that one may opt to average 
them over the structures /3, 

VMx) = (<(x)), VBBix) = {VU^)), VMx) = {vU^)). (3) 

Now that the pair potentials are known, the multi-body, non-pair ECI can be obtained 
with the conventional CEM procedure performed on the unrelaxed lattice, provided that all 
the pair potential contributions to the cohesive energy are subtracted out. 



Vaia] 



= E(^'')'^^'(«)' (4) 

where {C,~^)a is the (pseudo-) inverse of the correlation function matrix involving non-pair 
clusters a and structure /3. A£'^(a) the cohesive energy at a lattice parameter of a of the 
underlying lattice with the pair potential contributions subtracted out. A£'^(a) is derived 
from the right hand side of Eq. ([1]) , 

+ 2V:4B(R,,-K(l-yj], (5) 

where the occupation number p\ takes a value 1 if site i is occupied by the species indicated 
in the subscript and takes a value otherwise, which relates to the conventional Ising spin- 
like variable a as cij = 2p^ — 1 = 1 — 2p^. An interesting aspect of this hybrid cluster 
expansion is that unrelaxed cohesive energies are needed only - potentially a very significant 
saving in computational effort. 

Scheme II: Here, the mathematical form of the pair potential is known a priori, either 
due to the nature of the bonds or as a matter of expediency. The AA, BB, and AB pair 



potentials are determined through optimization of the adjustable parameters in the pair po- 
tential formula. For simple potentials especially, such as the Lennard- Jones potential,— these 
parameters are found quite simply requiring only a few calculations for unrelaxed structures 
of the pure elements and unrelaxed ordered structures of intermediate composition, or they 
could be determined from experimental data. The multi-body ECI are computedly lastly, 
in the same manner as under scheme I. 

Initial calculations with the first scheme at zero Kelvin have shown that NiAl with the 
unstable LIq structure properly relaxed to the B2 structure with a corresponding energy 
reduction of about 0.17 eV per atom and the ratio of the fee lattice parameters a/c changed 
from 1 to 1.414 as expected. In another preliminary test involving CuAu with initially the 
fcc-based LIq structure, the calculated lattice relaxation resulted in an energy decrease of 
0.03 eV/atom while the a/c ratio became 1.12, to be compared with the experimentally 
measured value of 1.07.— 

In conclusion, a hybrid pair potential - cluster expansion method has been developed 
which allows an efficient coupling of displacive and substitutional degrees of freedom in al- 
loys. The main merit of this hybrid cluster expansion is that it allows lattice relaxation to be 
modeled with relatively minor computational effort. It provides a more realistic energy func- 
tional for use in Monte Carlo, CVM and other lattice gas type simulations for thermodynamic 
properties of solids at finite temperatures. Unlike the conventional CEM which employs di- 
rectly ab initio relaxed energies, the hybrid CEM retains the degrees of freedom associated 
with relaxation in an explicit form so that relaxations can be extracted from lattice gas sim- 
ulations which use the hybrid energy functional. This means that temperature-dependent 
relaxations, and processes which are sensitive to local relaxations such as diffusion can be 
treated more realistically. Somewhat counter-intuitively, as the hybrid energy functional 
relies on non-relaxed structural energies it is computationally less demanding to derive than 
the conventional CEM. 
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